HD-THEP 98-40 



Time evolution of observables out of thermal equilibrium. 

Herbert Nachbagauer 
Institut fur Theoretische Physik, Universitat Heidelberg, 
Philosophenweg 16, D-69120 Heidelberg, Germany 
(Talk given at the 5th International Workshop on Thermal Field Theories 
and Their Applications, Regensburg, 10-14 August 1998) 



00 ■ 
On ■ 
On : 

> : 
o . 

£: 

\o : 

<N . 

(N ' 
> 

\Q . 
(N ■ 
On i 
O . 
00 ■ 

5 : 
^ : 

a: 
i 

<D . 
^ ■ 



X 



We propose a new approximation-technique to deal with 
the exact macroscopic integro-differential evolution equations 
of statistical systems which self-consistently accounts for dis- 
sipative effects. Concentrating on one and two point equal- 
time correlators, we develop the self-consistent method and 
apply it to a scalar field theory with quartic self-interaction. 
We derive the effective equations of motion and the corre- 
sponding macroscopic effective Hamiltonian. Non-locality in 
time appears in a natural way necessary to account for en- 
tropy generating processes. 



I. INTRODUCTION 

The question of how and if a statistical infinite di- 
mensional quantum system thermalises is one of the fun- 
damental problems of statistical mechanics, and has re- 
cently attracted much attention in connection with prob- 
lems in cosmology and relativistic heavy ion physics. We 
discussed some related issues in a recent work 0] for a 
zero dimensional quantum mechanical model. There, we 
concluded that a safe and consistent method to investi- 
gate on time evolution of macroscopic observables is to 
apply the projection operator technique. In this work, we 
systematically apply the new self-consistent approxima- 
tion technique anticipated in ffl to a scalar field theory 
with quartic self-interaction. For systematic reasons, let 
us therefore briefly review the line of arguments leading 
to that concept. 

In statistical theory, a level of observation is defined 
by a set of Hermitian so-called relevant operators the 
expectation values of which are the macroscopic observ- 
ables of the ensemble under consideration. The dynamics 
is defined by the microscopic Hamiltonian which enters in 
the exact integro-differential equations of motion for the 
macroscopic observables, as derived by the projection- 
operator technique |2j. These equations solve the com- 
plete BBGKY-hierarchy in integral form and are closed in 
the observables under consideration. Correlators not ex- 
pressible in terms of the basic observables are eliminated 
by the projection procedure but leave their imprints im- 
plicitly in the defining equation of a projected evolution 
operator. The problem of resolving the infinite set of 
differential equations of the BBGKY-hierarchy, and sub- 
sequent reduction to the macroscopically relevant quan- 



tities, has thus been converted into approximating that 
operator in a self-consistent way. Only in some particular 
cases, as for effectively non-interacting systems defined 
by a Hamiltonian which contain polynomials of at most 
second order in position and momentum operators, the 
projected time evolution operator can be integrated with- 
out explicit dependence on projectors. A side-condition 
for exact integrability is that also the level of observation 
must be chosen to be at most of second order (quadratic 
level of observation). The physical reason for exact in- 
tegrability is dynamical closure, i.e. the fact that time 
evolution driven by effectively non-interacting Hamilto- 
nians does not lead out of the linear space spanned by the 
operators of a quadratic level of observation. Typical ex- 
amples of this type are time dependent mass-parameters 
or QED without quantized fermions, but including clas- 
sical external sources. 

For the interacting case, the expansion of the projected 
evolution operator into a coupling parameter labeling the 
interaction part is discussed extensively in the literature 
However, that expansion scheme turns out to be 
inconsistent with the projected dynamics at second or- 
der even when we introduce an effective quadratic micro- 
dynamical Hamiltonian and expand the projected time 
evolution operator around it. That was demonstrated 
explicitly in a zero-dimensional toy model recently [Q. 
The break-down of that kind of expansions becomes ap- 
parent if one integrates the system for times large com- 
pared to the inverse coupling of the interaction part. 
The physical reason for that may best be illustrated by 
considering the uncertainty product expressible in terms 
of the quadratic observables. Even for effectively non- 
interacting Hamiltonians with arbitrary time dependent 
coefficients, the uncertainty - which essentially is the vol- 
ume of phase space cells - is an exact constant of motion. 
However, the dissipative part of the equations of motion 
for the observables does not respect constant uncertainty, 
and evolves it to numerically large values at secular time 
scales. On the other hand, increasing uncertainty and in- 
creasing phase space volume of the corresponding macro- 
scopic dynamics is an essential feature of non-equilibrium 
statistical systems. Moreover, the uncertainty product is 
directly related to the relevant entropy of a quadratic 
level of observation. We concluded that any effective mi- 
croscopic Hamiltonian is intrinsically incompatible with 
time variant entropy and cannot be used to approximate 
the projected evolution operator. One necessarily has to 
generalize to an effective Liouvillian acting in the space 
of relevant operators and expand around its associated 
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so-called reduced time evolution operator. 

The task of the present investigation is to derive the 
effective equations of motion by approximating the pro- 
jected evolution operator by its reduced counterpart to 
the first non-trivial order. We also construct the cor- 
responding effective Hamiltonian which generates the 
macroscopic equations of motion. 

II. DEFINITION OF THE PROBLEM 

(i) A mixed state of a quantum system (configuration), 
both in zero dimensional quantum mechanics as well as 
in quantum field theory, is described by a density opera- 
tor, which, in order to allow a probability interpretation, 
must be a hermitian trace-class operator with positive 
eigenvalues, ft is an intrinsic feature of quantum sys- 
tems that the density matrix is fictive in the sense that 
only its diagonal elements correspond to physical proba- 
bilities, the other dependencies being phases which enter 
in expectation values via interference effects^]. Alterna- 
tively, one may characterize a configuration by the ex- 
pectation values of hermitian operators. A generic set 
of those representing a complete set of observables is 
given by the mutually orthogonal hermitian projectors 
constructed from the eigenvectors of the density matrix. 
A complete set of (not necessarily commuting) observ- 
ables contains all information about the density matrix 
such that any observable can be expressed in terms of 
the complete set. 

(ii) The system may be assigned a dynamical struc- 
ture. Motion is defined as a sequence of possible states 
having a constant expectation value of the Hamiltonian 
operator, the energy of the configuration. Quantum me- 
chanical time parameterizes those configurations which 
are assumed to have time-independent probabilities and 
phases for autonomous systems. The time evolution gen- 
erated by the micro-dynamical Hamiltonian can be ex- 
pressed by first order differential operator-equations in 
time for the density matrix (von Neumann equation) in 
the Schrodinger representation, 

ij t p{t) = [H,p{t)], (1) 

together with an initial condition p(0). Hermitian 
Hamilton-operators generate unitary time evolution 
which in turn is necessary to allow for a probability in- 
terpretation. Once the problem of time evolution of the 
mixed state is solved for a given initial configuration, ob- 
servables can be calculated as expectation values from 
the evolved density matrix. 



(iii) The statistical description of a system is based on a 
reduction procedure selecting an in general much smaller 
number of (macroscopic) observed quantities out of the 
complete set of observables. This subset defines the level 
of observation (description), characterized by the set of 
relevant operators E — {F^}. The process of reduction 
from the density operator to the set of the so-called rel- 
evant quantities in general involves loss of information. 
Here, we will concentrate on levels of observation that 
once chosen, will be kept fixed during time evolution. 
That constraint can be relaxed too if necessary Q . 

(iv) The reduction can be dynamically trivial if it com- 
mutes with time evolution. In that case of dynamical 
closure, the Liouvillian maps the operators of the level 
of description on a linear combination of them. The ob- 
servables of a level of description at a certain time are 
sufficient to determine them at any time and evolution 
induces neither information gain nor information loss at 
the level of observation, the associated entropy being con- 
stant. On account of the canonical commutation relation, 
the most general Hamiltonian admitting a finite dimen- 
sional dynamically closed level of observation can contain 
constant, linear and quadratic expressions in position and 
momentum operators^]. The corresponding dynamically 
closed sets of observables correspond to sums of poly- 
nomials of finite order in the canonical operators. If the 
reduction is non-trivial, one can extend the level of obser- 
vation to render it trivial. That may involve an infinite 
number of operators in which case the system is a truly 
interacting one, and the only dynamically closed set of 
operators corresponds to a complete set of observables. 

(v) For truly interacting dynamical systems, the com- 
plete initial density matrix influences observables at later 
times. Its definition calls for an additional principle to 
construct it from the reduced set of initial data. Infor- 
mation theory proposes to apply Shannon's theory of 
entropy Q to that physical problem. Jaynes' principle 
of maximum entropy jq] fixes the generalized canonical 
density operator as initial condition. It can be shown 
not to contain more information than the initial set of 
observables. We want to point out that this choice is 
the statistically most probable, but the underlying con- 
cept of ensemble averages of identical systems evolving 
from variant initial preparations does not imply actual 
the preparation of the physical system in that state. 

The reduced description of a quantum system is 
achieved by passing from the density matrix to a statis- 
tical operator characteristic for the ensemble under con- 
sideration. The best guess for that operator compatible 
with the principle of maximum entropy has the functional 



2 Regardless of additional dependencies on other c-number 
quantities, including time and classical field strength, we call 
x That is the fundamental difference to classical phase space f ree Hamiltonian 

averages. 
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form 



ft = exp(— /i^Jv) 



It contains exactly the same amount of information on 
the system as the statistical observables g v = tr (pF v ) = 
tr (1ZJ- V ), and maximizes the relevant entropy functional 
of the level of observation, 



S £ = -tr (ft log ft) = p v g v 



(2) 



Reducing the description of a quantum mixed state to the 
evolution of the accompanying statistical operator ft re- 
quires to derive a modified evolution equation which does 
not lead out of the level of observation. The problem ba- 
sically boils down to introducing suitable projection op- 
erators and the study of their associated time evolution. 
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ternatively, one may solve the BBGKY-hierarchy which, 
however, is problematic since it involves an ad hoc trun- 
cation procedure. Problems with that approach have 
been discussed extensively elsewhere p| . 



III. PROJECTION OPERATOR TECHNIQUE 

We briefly review some basic terms of the projec- 
tion operator approach |||| serving as starting point 
of the present investigation. There, one defines the 
level of observation by the finite set of Hermitian op- 
erators £ — {F v } including J-q = TL for convenience. 
The accompanying canonical density operator defined 
by lZ(t) = exp(— p v (t)T v ) contains time dependent La- 
grange multipliers p u {t) which are functions of g v (t) such 
that g v (t) = ti{!F v p{t)) = tr (,F„ft.(t)). We emphasize 
that other than the density matrix pit), ft(i) in general 
does not evolve according to a Schrodinger equation. 

For the corresponding operator expectation values 
g v (t), the closed exact equation (Robertson equation) of 
motion is found to read 



-itr [T v LoTl{t)) 
dt 1 tr (T V L o T(t, t') o Q(f') o L o ft(t')) ■ (3) 



Here L stands for the Liouvillian induced by the micro- 
dynamical Hamiltonian, and the action of L on some op- 
erator X is defined by0 L o X = [H, X] . 

The r.h.s. of (|^) introduces two new symbols which call 
for further explanation. The projector Q(t) is defined by 



tx{OxQ{t) o Q 2 ) = tx{O x 2 ) 



dtr (OxU(t)) 
dg v {t) 



tr(JvO a ) 



(4) 



for any pair of operators 0i,02- In the particular case 
of 0\ being in the linear space of the level of observation 
£ , the rhs of Eq. (||) vanishes. 

The non-unitary projected evolution operator T(f , t') 
is a solution of 

d 

—T(t,t') = -iQ(t)oLoT(t,t') and 

^T(t,t') = iT(t,t')oQ(t')oL 

with initial condition T(t,t) = 1. It will turn out to 
be the key problem to construct this operator from its 



FIG. 1. Overview of relations on time evolution of mixed 
states and statistical time evolution. The arrows indicate the 
sense of deduction. 

An overview of the systematics is given in Fig. (1). 
Here, we adopt the Robertson equation approach. Al- 



3 The expression A o B in general denotes the action of op- 
erator A on operator B. There may be no representation of 
A o in terms of a commutator. 
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definition for the particular Hamiltonian under consider- 
ation. Then, all trace expressions at the r.h.s. of @ can 
at least in principle be expressed in terms of the <?„, and 
the system is a closed intcgro-diffcrential equation in the 
c- number expectation values. 

In this investigation, we will complete £ such that 
the action of the Hamiltonian induced Liouvillian can 
be split into L = Lq + Zj, and the free part be dy- 
namically closed with respect to the level of observation, 
LqoFv = [HqjJ-'v] = flv^J-^. In that case, in the integral 
part of Eq. (||), the complete Liouvillian can be replaced 
by Li. 

Expectation values of operators which cannot be rep- 
resented as linear combinations of the level of observation 
get additional contributions to their 7?.-averages, 

tx{o P (t)) = tx{on{t)) 

-i [ d«'tr(OT(M')oQ(i') i ^(i'))> (5) 
Jo 

where again the dynamically closed part in L does not 
contribute to the integral. 

IV. PROJECTED EVOLUTION 

The key problem is to find a consistent approximation 
scheme to solve the defining equation for the projected 
evolution operator T(t,t'). That scheme — order by or- 
der - must (i) be compatible with the Robertson equation 
and (ii) conserve exact constants of motion. In particu- 
lar, the expectation value of the Hamiltonian operator 
being a generic constant of motion must be conserved 
which is a non-trivial condition since, due to relation 
, that generally involves integral terms non-local in 
time. The query for a consistent approximation scheme 
naturally raises the question whether there exists an ef- 
fective micro-canonical Hamiltonian which may serve as 
a zeroth-order approximation being a staring point of an 
expansion. 

A fist approach to that problem is to study the small A 
expansion of H = Hq + XHj . If one evaluates the corre- 
sponding expressions to second order in the coupling and 
solves the resulting system of integro-differential equa- 
tions of motion, one finds that the expansion breaks down 
at typical secular time scales of the order t ~ 1/A pj. 

An improvement of that expansion scheme in equilib- 
rium theory can be achieved by resummation. In our 
case, that amounts to introduce an effective Hamiltonian. 
In particular, one may try to make an quadratic Ansatz 
of the form 

H eS = a(t; x, y)U(x)U(y) + b(t; x, y)${x)$(y) + ... 

and determine the time dependent coefficients by self- 
consistency requirements. However, it turns out that 
even that improved scheme fails to meet the consistency 
conditions (i,ii) at secular time scales Q. The physical 



reason for that is the presence of memory effects. In 
particular, although the effective Hamiltonian correctly 
describes time evolution for short times t <C 1 / A in an ap- 
propriate way, the time evolved quantity H e Q does not 
coincide with effective Hamiltonian constructed by the 
consistency requirements at later times. We conclude 
that only the enveloping dynamics is an effective Hamil- 
tonian one. Furthermore, effective Hamiltonian dynam- 
ics automatically gives rise to an additional constant of 
motion related to the uncertainty product of the config- 
uration. On account of the uncertainty-entropy relation 
for quadratic observables, that results in conservation of 
the relevant entropy contradicting constant information 
loss inherent in the dissipative time evolution of statis- 
tical systems. Also, locality in time conflicts with the 
presence of memory terms, which get important at large 
time scales. Third, an effective Hamiltonian automati- 
cally induces a unitary effective time evolution operator 
which again conserves observable information. The con- 
clusion is that there generally exists no effective micro- 
dynamical Hamiltonian dynamics which is compatible 
with the Robertson-equation and conservation of energy. 
We thus have to completely abandon the query for micro- 
dynamical Hamiltonians in favor of a the more general 
concept of Liouvillian dynamics. 

The key idea is to formally introduce an effective evolu- 
tion operator which evolves the accompanying statistical 
operator lZ(t) in a form-invariant way and expand the 
projected evolution operator around that reduced time 
evolution. 

The reduced time evolution super-operator V(t, t') is 
uniquely defined by two demands. Firstly, it does not 
lead out of the linear space of the relevant operators (ob- 
servables). Secondly, it evolves the statistical operator, 

Tl{t) = V(t,t') oTZ(t'). (6) 

That condition already is very restrictive. Due to the 
particular form of the accompanying statistical operator, 
the differentiated equation suggests to write a solution as 
a time-ordered exponential 

V(t,t')=T exp (-i £ drL(r) o V (7) 

Since the derivative d t obeys the Leibnitz product-rule, 
the effective Liouvillian-Operator L{t) has to have the 
algebraic property of a derivation, L(t) o (AB) = (L(t) o 
A)B + A(L(t) o B). We want to point out that in general 
this Liouvillian is not an anti-hermitian operator, L =/= 
—L > and thus cannot be cast as a commutator with some 
effective Hamiltonian, L ^ [H, .]. 

Furthermore, on account of linearity in the level of ob- 
servation, the reduced time evolution can be expressed 
in terms of a matrix evolution operator, V(t,f) o T a — 
V a f3(t,t')J-f3. The reduced time evolution operator still 
possesses the property of transitivity, V(t, t') o V(t' , t") = 
V(t,t")- 
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In practical calculations, it turns out to be convenient 
to act with the reduced time evolution operator to its 
left which amounts to define the adjoint action induced 
by the trace as inner product by 

tr (OV(t, t') o 1l{t')) = tr ((V f (t, t') o 0)K(t')) 

= tr{0{t,t')TZ(t')). (8) 

The adjoint evolved operator 0(t,t') is the general- 
ization of the free Heisenberg operator to O for non- 
Hamiltonian time evolution. Due to the particular prop- 
erties of V(t, i'), one can also express time evolution of the 
generalized Heisenberg operators of the basis by means 
of a matrix multiplication, .F Q (i,i') = V^(t, t') o T a — 
V a g{t, t'^Fp. The matrix elements are restricted by the 
conditions 

9a(t)=Vt f3 (t,t')g l3 (t'). (9) 
In addition to that, the property 

V f (t, t') o {T a T ) = (V* (t, t')oT a ) (V f (t, t') o T p ) (10) 

which is a direct consequence of L{t) being a deriva- 
tion, further restricts the number of independent ele- 
ments in the evolution matrix if some symmetrised prod- 
ucts {!F a ,!Ffj} are also in the operator basis of the level 
of observation. In general, however, those conditions are 
not sufficient to fix the evolution matrix completely. 

The self-consistent method now requires to expand the 
projected evolution operator in terms of the reduced time 
evolution operator. It is part of our strategy to exploit 
the remaining freedom in the evolution matrix to satisfy 
consistency with the Robertson evolution equations, and 
possible conserved quantities. As will be shown in the 
case of the local and bilocal equal-time correlators as 
observables, that fixes the dissipative evolution matrix 
in a consistent way. 

V. REDUCED EVOLUTION OPERATOR FOR 
ONE AND TWO POINT CORRELATORS 

In order to define the ensemble, a particular set of 
observables must be chosen. Among various choices of 
the level of observation, the special case of quadratic ob- 
servables plays a particular role. This set of Hermitian 
operators that can be constructed out of the field op- 
erators and its canonically conjugate momentum, {1 , 
$(x), n(x), Q(x)$(y), W(x, y) = *(x)IL(y) + Il(y)$(x), 
U(x)H(y)}, and give rise to a quasi Gaussian bilocal ac- 
companying statistical operator. This particular choice 
of the level of observation renders possible to apply a gen- 
eralized Wick theorem in the evaluation of weightened 
quantities (0) t — tr (01Z(i)). Then expectation values 
of polynomials in <E>,n can always be expressed as prod- 
ucts of contracted pairs of the macroscopic observables 
<p(x\t) — (<&(x)) t , n(x\t) — (II(x)) t , and the correlators 



<p 2 (x, y\t) = (*(a:)*(i/)) t , tt 2 (x, y\t) = {U{x)n(y)) t , 
w(x,y\t) = (W(x,y)) t . (11) 

We now construct the reduced time evolution operator 
for this level of observation. Since 1 , $, II already rep- 
resent a basis out of which the quadratic quantities are 
formed, it suffices to investigate on the evolution matrix 
elements contained in 

( SSlo ) = ^c*,o - ( S ) +-yc*,0 (12) 

Instead of considering the three vector of operators 
($,II, 1 ), we have eliminated the trivial equation for 
the unity which leaves us with an inhomogeneity. The 
elements of K(t, t') and -f(t, t') contain the remaining in- 
dependent elements of V'(i, i'). The coefficients of the 
matrix K and the vector 7 are constraint by (||), to wit, 

)-*<«.<>■(!$} < i3 > 

For briefness, we will use the symbolic matrix nota- 
tions {A-B){x,z) = f y A(x,y)B(y,z), (A ■ v)(x) = 

J y A ( x >y) v (y), (A T )(x,y) = A(y,x) and (v®w)(x,y) = 
v(x)w(y) whenever the meaning is unambiguous. If we 
shift the operators <&, II by their expectation values, the 
constant 7 can be eliminated, and the product condition 
( |To| ) gives rise to the relation 

Z{t) = K(t, t') ■ Z(t') ■ K T (t, t'), (14) 

where Z contains the quantum widths, 

z ._ ( <p2-<p®<p ±w-p®Tr\ (15) 
\_ — tt ® (p ir 2 — it ®ir J 

In order to solve the condition (fl4|), we remember that 
the reduced time evolution operator V(t, t') is a time or- 
dered product which by construction obeys a transitivity 
condition V (*,*') = V(t, t")V(t", t'). If we let t" = it 
becomes apparent that double time dependence factor- 
izes into a product where the second factor is the inverse 
first operator at t' . Consequently, we can also split the 
Heisenberg coefficients into 

K(t,t') = A(t)- A- 1 ^'). (16) 

If we put back that form into the restriction (|l4|), one 
finds that the quantum widths factorize into 

Z{t) = A{t)-A T {t). (17) 

We emphasize that Eq. jl7| ) does not put any restriction 
on the free choice of the initial expectation values. Z by 
construction being a symmetric matrix, can always be 
diagonalized by a unitary transformation and thus can 
always be written asZ = RDR T = Ry/D(Ry/D) T . 

Further conditions can be obtained from compatibility 
of the time evolution ([l3]) with the equations of motion. 
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For the very general class of Hamiltonians composed of a 
quadratic kinetic part II 2 /2 and a canonical momentum 
independent potential V($), the equations of motion (||) 
for the position coordinate are exactly <fi(t) = n(t) and 
do not get any further dissipative corrections. It follows 
in this case that A has to have the structure 



containing the positive definite local operator O z = m 2 - 
A z , and interaction 



A{t) 



a(t) 0(t) 
a(t) p(t) 



(18) 



Together with ( p~T|) that form also implies 2(f2 (t) = 
w(t) + w T {t) which is in fact also exactly valid for theo- 
ries with quadratic kinetic terms. For Hamiltonians with 
a different kinetic term, one has to keep a more general 
form of the matrix A(t) which renders the algebra more 
complicated, but the construction still works along the 
same lines. 

A solution of the Eqs. (p7|-^8[) can be parameterized 
by a(t) + i(3(t) =X(t) ■ U(t), where U(t) has to be a uni- 
tary matrix and X(t) be real. When plugging in that, 
we benefit from the fact that unitary matrices can be re- 
garded as solutions of the first order differential equation 
U(t) = ifl(t) ■ U(t) with n(t) = fl T (t) being real. We get 

(f 2 - <p® f = X ■ X T , -w - ip <E> 7T = X ■ P T (19) 

with P = X, and 

n 2 -Tr®Tr-P-P T = X-n-rL-X T . (20) 

These relations can be combined into the matrix evolu- 
tion operator, which first row elements are found to read 



a{t,t') =X(t)-C(t,t')-X-\t') - 
0{t,t')-P{t')-X-\t') 



(21) 



and 



f3(t,t') = X{t) ■ S(t,t') ■ n^it') ■ X-\t'). (22) 



The functional real generalisation C,S of the trigono- 
metric functions emerge from biquadratic expressions of 
the real and imaginary part of the matrix U(t) and solve 
the differential equations 

C(t,t') = -Q(t) -S(t,t'), S(t,t') = n(t)-C(t,t') (23) 

with initial conditions C(t,t) = 1, S(t,t) = 0. The 
complete form of the matrix elements of the reduced evo- 
lution operator can now be expressed in terms of the 
quadratic observables and one action variable f2. 



VI. EQUATIONS OF MOTION 

To be specific, let us study a model with Hamiltonian 
operator H = Ho + Hj with dynamically free part 

i/ = J|(n(z) 2 +$(z)O 2 <D(z)) (24) 



Hi = \l <!>V... 



(25) 



For that particular Hamiltonian, the equations of motion 
(H) evaluate to the macroscopic evolution equations 

V(x)=n(x), (26) 
•?t(x) = —<> x Lp(x) + 2\(p(x) (2ip 2 (x) — 3(f2(x, x)j 

+A 2 E 7r (x), (27) 

•i>2(x,y) = -(w(x,y) + w(y,x)), (28) 

ttsOe, V) = -7}<>xw(x, y) + A (4ip 3 (x)ir(y)- 

3(f 2 (x, x)w(x, y)) + (x <-> y) + 
X 2 ^ 2 (x,y), (29) 
w(x,y) = 2ir 2 (x,y) - (O x + O y )ip 2 {x,y) 

+4\(2v(x)cp 3 {y)-3<p 2 {x,y)(p2(y,y)) 
+\ 2 ^ w {x 1 y). (30) 

where the expressions £„ contain the dissipative integral 
part. Analogously, the expectation value of the Hamilto- 
nian is found to read 



(H) t = 

i J (n 2 {z,z) + ^O z i P2 (z,z) + \(3^ 2 2 {z,z)~2i P i (z)) 



+A 2 S ff (z)) 



(31) 



The evaluation of the integral terms requires an ap- 
proximation of the projected evolution operator T(t,t'). 
By making the replacement, T{t,t') — > V(t,t') in we 
apply the first term of a systematic approximation of the 
complete Liouvillian L around the reduced evolution gen- 
erated by Lit). After acting with the conjugate of V on 
the operators to its left which transforms them into the 
generalized Heisenberg picture, we are left with averaged 
operator products which can be split into pairs by virtue 
of the generalized Wick theorem. There, the propagator 

G(x, y\t, t') = 2(X(t) ■ C(t, t') ■ X T (t'))(x, y) 
= ([$(x\t,t')-<p(x\t)][<Z>(y)-<p(y\t')] + h.c.) tl (32) 

appears in a natural way. 

The dissipative expressions have the form of memory 
integrals, 



E„(i)= / / dt'a u {t,t') (33) 

z 

where, after some lengthy algebra, the integrands are 
found to read 



G 



(34) 



a w (x) = 6/3(1, (z) (3G 2 {x, z) - (3 2 (x, z)) 

3 \$(y, z)G{x, z) (G 2 (x, z) - 3/3 2 (ir, z)) + 
^(x, z)G(y, z) (3G 2 (a;, z) - /? 2 (ie, z)) + 
6(/3(y,z)(G 2 ( 2 ;,z)-/3 2 ( a ;,z)) + 

2f3{x,z)G{x,z)G{y,zj) ip(x)ip'(z)\ + 
a w (x)ip(y) + (x <-> y), (35) 

6 [jS(y, z)G(x, z) (3G 2 (y, z) - /3 2 (y, z)) + 
0(x,z)G( y> z) (G 2 (y,z)-3[3 2 (y,z)) + 
2/3(1/, *) (3G 2 (y, z) - /3 2 (y, z)) ^(x)^'(z) + 
6 (/?(*, z)G 2 (y, z) + 2/%, z)G(x, z)G(y, z) 

- / 9(x,z)/3 2 ( 2 /,z)) + 
6 (/%, z)G(x, z) - /3(x, z)G(y, z)) cp 2 (y, y)ip' 2 (z, z)] , 

(36) 

<T H (x) = 

-3/?(ar, z) (G(z, z)(G 2 (x, z) - (3 2 (x, z))+ 
2(,5(x)^'(z)(3G 2 (z, a) - /3 2 (x, z))) (37) 

with ip'(x) = ip(x\t'), (p 2 (x,y) = tp(x,y\t') and other 
dependencies on time arguments have been suppressed. 
The Eqs. ( p(|]3"o| ) together with these relations finally 
form the closed set of macroscopic equations of motion. 

VII. EFFECTIVE HAMILTONIAN 

A natural question to ask is if the effective equations 
of motion which contain a memory term non-local in 
time still can be generated by a macroscopic Hamiltonian 
or some action principle. A straightforward candidate 
for that Hamiltonian is the expectation value h = (H) t 
which, however, has to be recast in terms of canonically 
conjugate variables. 

Leaving aside the dissipative contributions for the mo- 
ment, the effective Hamiltonian was constructed system- 
atically recently |Q for the zero-dimensional case. In 
analogy to that result, the equation of motion ( p6|) sug- 
gest to introduce the macroscopically canonically con- 
jugate pair {tp(t), 7r(t)}. Furthermore, if we chose the 
positional uncertainty X(x,y\t) as bilocal second posi- 
tion coordinate, the corresponding canonically conjugate 
pair is found to be {X(x,y\t), P(x,y\t)}. Together with 
the relations ( ^9|j2C| ), we find that the non-dissipative 
part of the energy is composed of three contributions, 
hi = T + Vp t. + Vam.i with a kinetic term 

T=\j {7r(z)n(z) + (P-P T )(z,z)), (38) 

z 

a potential 



Vpot. — 



y(z)O z y(z) + -O z (X ■ X T ){z, z)+ 



A (3(A • X T ){z, zf + 6<p 2 {z)(X ■ X T ){z, z)+ 
V\z))) , (39) 

and an angular momentum term (A = X T ■ X • O), 



Km. =\J ((X T 1 ■ A - A T ■ A- X )(Z, Z)) . (40) 



That term comes from the kinetic part of (H) t when 
transforming into canonical coordinates jl],^) 

The Hamiltonian equations of motion for the variables 
ip, 7r; X, P as being generated by h\ are equivalent to the 
non-dissipative parts of the original Eqs. (^ , ^o| ). This 
system, however, has the additional integral of motion 
A. In order to reproduce that property by the effective 
Hamiltonian too, we introduce a canonically conjugate 
variable for the momentum A. The equations of motion 
for the conjugate pair {*f?(t), A(t)} then in fact yields a 
constant uncertainty A(t) = since the non-dissipative 
Hamiltonian is independent of 

This canonical structure can be carried over to the dis- 
sipative corrections. To see that in detail, we first realize 
that the energy corrections do not depend on the 
momentum A(t), but only on the quantities ip, 7r, X, P, S 
and C. Regarding the expectation value for the energy 
as functional of those canonical variables, it follows that 



SA 



(41) 



By comparison with the definition (|23j), it is possible to 
split off the factor if we regard the generalized trigono- 
metric functions as functionals of \& being subject to the 
defining equations 



6C[V](u,v) 

5^(x, y) 
5S[V](u,v) 



= -S(x,u)S(y,v), 
= -S(x,u)C(y,v). 



(42) 



with initial conditions C[* = 0] = 5, <S[* = 0] = 0. The 
complete macroscopic Hamiltonian equations of motion 
with canonical pairs {ip, 7r; X, P; 'J, .4} are thus generated 
by the classical effective macroscopic energy functional 



h[<p,w;X,P;*,A}=T+V ar , 



pot. 



A 2 E 



H 



(43) 



and are in fact equivalent to the original equations of mo- 
tion @Q. The effective macroscopic Hamilton func- 
tional is non-local in time, and accounts for the full his- 
tory of the system. The corresponding dynamical initial 
value problem is only defined for the initial time t = 
where the initial preparation is determined by Janynes' 
principle of maximum uncertainty. The Hamiltonian h, 
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however, cannot be generated as classical limit of a cor- 
responding microdynamical Hamiltonian, since the re- 
duced microdynamics is a Liouvillian one. We further 
remark that the relevant entropy defined in (|^) is auto- 
matically larger or equal at later times than its value at 
t = 0, though intermediate entropy fluctuations are not 
excluded from dynamics. We mention that by means of a 
standard Legendre transform one straightforwardly gets 
the corresponding macroscopic Lagrangian and effective 
action. 



VIII. CONCLUSION 

We constructed the effective macroscopic equations of 
motion and Hamiltonian for the statistical ensemble de- 
fined by expectation values and equal time correlators 
of position and momentum operators in a scalar theory 
with quartic interaction potential. The systematic con- 
struction based on the exact reduced equations of mo- 
tion for the corresponding level of observation have to 
be approximated self-consistently in order to avoid viola- 
tion of energy conservation and for compatibility with the 
Robertson-equation of motion. We solved that problem 
by a systematic approximation method based on an effec- 
tive microdynamic Liouvillian evolution compatible with 
time evolution of the accompanying statistical operator, 
rather than using an effective microdynamical Hamilto- 
nian which fails to satisfy the consistency requirements. 
The method fully accounts for entropy generation and 
thus can be regarded as suitable tool to study the ef- 
fects of dissipative time evolution even at arbitrarily large 
time scales. Moreover, the method intrinsically is supe- 
rior to the study of the BBGKY-hierarchy of equations 
of motion since it both includes a systematic principle 
to solve the initial preparation problem, and avoids the 
unresolved difficulty of how to truncate the BBGKY- 
hierarchy without implicitly introducing approximation 
method artifacts. Further applicaltions are planned for 
future research. 



[2] J. Rau and B. Miiller, Phys. Rep. 272 (1996) 1, and refer- 
ences therein; 

E. Fick and G. Sauermann, The Quantum Statistics of 
Dynamic Processes, Springer series in solid state sciences, 
86, (1990); 

K. H. Li, Phys. Rep. 134 (1986) 1; 

H. Grabert, Projection Operator Techniques in Statistical 
Mechanics, (Springer, Berlin 1982). 
[3] R. Balian, Y. Alhassid and H. Reinhardt, Phys. Rep. 131 
(1986) 1. 

[4] C. E. Shannon and W. Weaver, The Mathematical Theory 
of Communication, (University of Illinois Press, Urbana, 
1949). 

[5] E. T. Jaynes, Papers on Probability, Statistics and Statis- 
tical Physics, (Kluwer Academic, Dordrecht, 1989), ed. by 
R. D. Rosenkrantz, in part. Phys. Rev. 106 (1957) 620. 

[6] F. Cooper, S. Habib, Y. Kluger and E. Mottola, 
Phys. Rev. D55 (1997) 6471. 



ACKNOWLEDGMENTS 



I would like to thank U. Heinz and the local organizers 
of the "5th International Workshop on Thermal Field 
Theories and Their Applications" for all their efforts. 



[1] H. Nachbagauer, Dissipative Time Evolution of Observ- 
ables in N on- equilibrium Statistical Quant um Systems, 
preprint, HD-THEP-98-32, |hep-th/9807144 To be pub- 
lished in Euro. Phys. J. C. 



8 



